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DETERMINATION OF PITCH AND ROLL 
FROM HORIZON EXPOSURES ON KH-4 SATELLITE PHOTOGRAPHY 


INTRODUCTION 


Horizon photographs are used to determine 
the pitch and roll components of the KH-4 
vehicle's attitude or, more properly, of an or- 
biting camera's attitude. Horizon photographs 
do not give any vaw information in cases where 
the primary or pan camera is aimed atthe nadir 
(Figure 1). The determination of pitch and roll 
is simple and straightforward, and good approxi- 
mations of both pitch and roll are easily obtained 
from simple measurements and computation. * 
When the primary camera is not oriented toward 
the nadir and the horizon cameras are attached 
to the primary camera sothat the horizon photos 
have an apparent swing, as is the case with KH-4, 
the measurements are restricted to a smaller 
section of the horizon image and the data re- 
duction becomes more complicated (Figure 2). 









FIGURE 1. HORIZON IMAGE FOR VERTICALLY ORIENTED 
PRIMARY CAMERA. 


Many reduction systems have been pro- 
posed and used. They range from the ex- 
tremely simple ones where simple scales are 
used to measure the image and the computation 
is done graphically, to the complex ones where 
precise machines are used to measure the image 
and the computation is accomplished with a 
computer. The following is presented as a 
solution falling between the extremes, but def- 
initely toward the more complex approach. 
This determination of pitch and roll from 
horizon photographs is tailored to be practical 
for the reduction of a large number of photo- 
graphs. it is assumed that measurements are 
made with a comparator, such as a Mann 
instrument, and that the computation is per- 
formed by a computer. 








FIGURE 2. 
CAMERA. 


HORIZON IMAGE FOR PITCHED PRIMARY 


DEFINITIONS 


Yaw, pitch, and roll are taken as orthogonal 
rotations about the three mutually perpendicular 
primary photo axes having their origin at the lens 
(Figure 3). Yaw is considered primary, pitch 
secondary, and roll tertiary. Since a mechanical 


*The final result is not error free: a significant error arises 
in the identification and measurement of the horizon image. 


gimbal system is not involved, the ordering of the 
orientation rotations is arbitrary though neces- 
sary to adequately define them. Swing, S, is a 
rotation of the x-y-axes about the z-axis in 
the plane of the photograph. Swing is positive 
when the x-y-axes are rotated tothe left as seen 
from a point on the positive z-axis. 
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FIGURE 3. ATTITUDE COMPONENTS. 


MEASUREMENTS 


The measurement of image points is done 
with respect to an axis system onthe photograph 
having its origin at the principal point. Ten 
points are measured on each horizon frame, as 
shown in Figures 4 and 5. The4 fiducial marks 
are read first and then 6 points on the horizon 
image. Points 5 and 10 should be chosen as 





FIGURE 4. MEASUREMENT POINTS. 


29:2 


Handle Vic 


Control Only 


we wit! 


far apart as possible, and points 6 through 9 
should be approximately evenly spaced between 
them. 


COMPUTATION 


Assume the following photo measurements: 


p, (x,,y,) = end fiducial mark 
Pa (xy, ¥,)= Side fiducial mark 
Ps (Xs, ¥3)= Opposite end fiducial mark 
p, (*,,y,)= Opposite side fiducial mark 
p, through p,;, = approximately evenly spaced 
points on the horizon image. 
The principal point is determined from the inter- 
section of two straight lines. Instrument toler-. 
ances in placing the fiducial marks prohibit the 
method of taking the average coordinates of op- 
posite fiducial marks. 
straight line: 
yo( ESB) BEE) 
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From the equation of a 








These equations represent the fiducial lines de- 
fined by the opposite fiducial marks. Their 
simultaneous solution gives the coordinates of 
the principal point. It is convenient atthis step, 
but not necessary, to reduce the coordinates of 
the points on the horizon curve to the principal 
point origin, i.e.: let 
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FIGURE 5. ROLL COMPONENT (d) AND SWING (S). 
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Assuming the horizon image to be circular, 
the center is found from the circle of best fit in 
the sense of least squares. It is worthwhile to 
note that the horizon image is not circular but, 
because such a small section of the actual hyper- 
bolic curve is given, there is no practical dif- 


ference. It is convenient to find the center of 
the circle. From the general equation of a 
circle: 


x74 y2+ Ax + BY+C=O0. 


There are six observation equations of the form 
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which give the following normal equations: 
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and 


x y 1 1 
+ |——| + [| ————]} + =O. 
(x? + y3)? (x? + y4)? (x? + y 9)? x2+y? 


The simultaneous solution of these normal equa- 
tions gives the coefficients of the equation of the 
circle. Then the center coordinates of the circle 
are: 























x =# and Ye = 8. 
The roll component, 4, onthe photo is merely the 
difference between the radius of the circle and 


the distance fromthe principal point to the center 
of the circle (Figure 5): 


d=h (Vat ~ VaA?+B2-4C 


The roll component p , for the horizon camera is 
determined bv: 


d 
stan) =, 
c { 


where | is the calibrated focal length of the 
horizon camera. 
Swing is given by: 
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S' = ta 


Because the y-axes of the horizon camera 
and the primarv camera are not exactly parallel, 


the calibration data are used to rotate the hori- 
zon camera axes to parallel those of the primary 
camera. This rotation is significant in swing. 
The correction to the swing is 
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D 
where pb, N, and k arecalibration numbers given 


in the KH~4 calibration data for the port horizon 
Then 


AS=t 


camera of the forward set. 
S=S'+ AS. 
The computation thus far is computed for 
both the starboard and port horizon cameras 
and then the average roll component, 
| a ee 
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is used to compute the dipangle, p',for the port 
horizon camera measured inthe principal plane. 
This operation begins an iterative loop which 
usually yields good results in three cycles. The 
projection of the preset or "built-in" dip, p, to 
the principal plane is 

D'=ain™ (sin Dcoa gd), .....-. (1) 
where ¢ is replaced by S for the first approxi- 


mation. Then the roll in the principal plane is 
Pia D'+ Pm, ......26.4. (2) 
and the roll in the pitched plane is 
—_ sin?’ 
Pena! (Se | i agttes fon ans (3) 


where ¢ is replaced bv S for the first approxi- 
mation. The azimuth of the principal plane is 
computed from the expression 

behing ara a €4) 
where ¢ is replaced by S for the first approxi- 
mation. Delta is really half the tilt of the 


horizon photo, therefore 


« =tan~! (tan P ain d), 


Sa a owed ere (oD 
which gives all of the components required to 
compute the pitch: 

& =tan7! [(tan S + tan «) cos = tan 5]. ae (6) 


New values for °, 6, «, and 4 are computed 
using the last computed values in equations 
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(1) through (6). Three iterations are usually 
enough. The values of (equation3) and (equa- 
tion 6) refer to the left horizon photo. 

The pitch and roll for the primary camera is 
determined bv adding the calibration data. The 
KH-4 calibration data provides the following: 

1. yy. = the y-coordinate of the pan camera 
Then the pitch for the pan 


i 


calibration mark. 
camera is given by 


where F is the calibrated focal length of the 
pan camera. In the case of roll, the KH-4 
calibration data is: 
2. \, = the x-coordinate of the horizon cam- 
era calibration mark, which is co-planar with 
3. X,, = the x-coordinate of the pan camera 


to the corresponding mark. 
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FIGURE 6. GEOMETRY OF TILTED HORIZON PHOTO. 
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Then the roll for the pan camera is 
Py =P + 73-90 + tan! At + cE 
: ane 
=P 15+ tan! at - ee Bn Siatta be Sg cine (8) 
If only one horizon photo is used, the al- 
titude is used in the computation. In this case 
e'’ is computedas above, and then tilt is deter- 


mined from R 
T=sin"! EE AG (9) 
R-H 





where R is the mean radius of the earth and H 
is the altitude above the surface. The insen- 
sitivity of this solution indicates that the radius 
approximation is valid. Then 


and 
Now the three-cycle iteration with equations 


(1) through (6) will vield the pitch and roll, for 
the horizon camera. 
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APPENDIX 


Derivation of Equations 3, 4, and 6 (Fig- 3. Equation (6): 


ure 6). 
) Gk _KQ-GQ 
1, Equation (3): SOR a ape? 
, P'B fsinPcos KQ = KM sec «=f tang sec «, 


sinP' =——= , 
f f GQ = GN tan « =f tan dtane, 








(225) & ftan @ sec = -f tan & tan « 
P =sin- tan § = 
sin cos d . Ttan 5 
2. Equation (4): Rearranging: 
p'p" fsinP sind tad fan Wane ten S + tan 6x80 « 
tan « = ee ENTS RARE NRE TAR e 
‘ P''C f cos P sec « 
« =tan™! (tan P sin g). @ = tan"'{an S + tan «) Cos « tan 5]. 
e 
= 5:5 
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